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ABSTRACT 

A planet of low mass orbiting in a two-dimensional gaseous disc generates a one-armed 
spiral wake. We explain this phenomenon as the result of constructive interference 
between wave modes in the disc, somewhat similar to the Kelvin wedge produced in 
the wake of a ship. The same feature is not expected in a three-dimensional disc with 
thermal stratification. 

Key words: accretion, accretion discs - hydrodynamics ~ planets and satellites: 
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1 INTRODUCTION 

The dynamical interaction between a planet and the disc 
in which it forms has important consequences for the evo- 
lution of the orbital elements of the planet. In particular, 
any imbalance between the torques exerted on the planet 
by the parts of the disc interior and exterior to its orbit 
causes the planet to migrate radially through the disc (e.g. 
Ward 1986). The planet-disc interaction has been studied 
using linear theory and, more recently, non-linear numerical 
simulations. Most of this work uses a two-dimensional de- 
scription that ignores the vertical structure and motion of 
the disc. 

In the linear theory (Goldreich & Tremaine 1979) the 
planetary potential is analysed into Fourier components, 
each having a definite azimuthal wavenumber and angu- 
lar pattern speed. For each Fourier mode, the response of 
the disc is determined by an inhomogeneous linear wave 
equation that includes resonances at various radii. Angu- 
lar momentum is exchanged between the planet and the 
disc mainly in the vicinity of Lindblad resonances (at which 
waves are launched) and the corotation resonance (where 
angular momentum is deposited). Goldreich & Tremaine 
(1979) provided analytical formulae for these torques, while 
Korycansky & Pollack (1993) determined them by solving 
the linear wave equations numerically. 

Neither of these studies was concerned with the form 
of the disturbance generated in the disc, which is obscured 
in the Fourier representation. In contrast, numerical simula- 
tions of the planet-disc interaction (e.g. Bryden et al. 1999; 
Kley 1999; Lubow, Seibert & Artymowicz 1999) aim to solve 
the non-linear fluid dynamical equations in real space. Most 
such simulations have treated the case of a Jovian planet, 
for which the planet-disc interaction is highly non-linear, 
especially close to the planet where a gap is opened in the 
disc. Further away, spiral shock waves are distinctly seen. 



A detailed comparison between linear and non-linear 
theory has been made by Miyoshi et al. (1999), who com- 
pared the torques generated in a local, three-dimensional 
simulation (an isothermal, shearing-sheet model) with those 
derived by solving the linear wave equations. They did not 
compare the form of the disturbance in real space, however. 

For planets of lower mass, a gap is not opened and the 
response of the disc should be well described by linear the- 
ory. Simulations for a terrestrial planet (Artymowicz 2001) 
indicate that a remarkable one-armed spiral wake is formed, 
but there is no evidence for features associated with indi- 
vidual Lindblad resonances. The relation between the wake 
and the linear Fourier analysis has been obscure. In this 
letter we explain this relation by showing how the wake is 
formed from the superposition of individual Fourier modes. 
We also discuss the differences that are to be expected in 
future three-dimensional simulations. 



2 LINEAR HYDRODYNAMIC WAVES IN A 
TWO-DIMENSIONAL DISC 

2.1 Dispersion relation and pitch angle 

Consider a two-dimensional gaseous disc around a star of 
mass M. Let the disc have angular velocity n(r), surface 
density E(r) and adiabatic sound speed Vs(r). Suppose that 
a planet of mass Mp orbits within the disc on a circular 
orbit of radius Vp and angular velocity flp . The planet inter- 
acts gravitationally with the disc and excites hydrodynamic 
waves having angular pattern speed Qp. We assume that the 
waves may be described using linear theory. 

For linear waves, any wave quantity X can be written 
in the form 



X(r,(t),t) = Re{X(r) exp [i<l>„(r-, i)] } , 



(1) 



where X is an amplitude that varies slowly with r, while 
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fc(r) dr + m(<f> - Qpt) 



(2) 



is a phase that varies rapidly with r. Here k is the radial 
wavenumber, which is real in regions of space where the 
wave propagates, and m is the azimuthal wavenumber. We 
consider only non-axisymmetric waves, and take m to be a 
positive integer by convention. We initially leave the con- 
stant of integration in equation (^) indefinite. 

The dispersion relation for tightly wound hydrodynamic 
waves in a two-dimensional gaseous disc is 



(3) 



where n{r) is the epicyclic frequency. We neglect the effects 
of self-gravitation and dissipative processes. The wavefronts, 
or lines of constant phase, are spirals defined by the equation 

^ = _A. (4) 
dr m 

We note that, if the term in equation may be neglected 
compared to the other terms, the pitch angle arctan(m/fcr) 
is independent of m. This raises the possibility, in principle, 
that constructive interference may occur between waves of 
different azimuthal wavenumber. Since k and Q, are usually 
of comparable magnitude (indeed, they are equal for a Kep- 
lerian disc), the term is relatively small when m is large, 
except close to the corotation radius where = f2p. 

2.2 Disc model 

We assume that the disc is Keplerian, so thatn 



K — Q, — 



The corotation radius is then 

1/3 



We assume further that the sound speed is given by 



(5) 



(6) 



(7) 



where e is a constant. This convenient assumption, which is 
not critical to the analysis that follows, implies that the disc 
has a constant angular semi-thickness H/r proportional to 
e, and a constant Mach number e~^. 

Adopting units such that GM = Sip = Tc = 1, we find 



e 

where 

r± -- 



( 3/2 3/2\ / 



^3/2 _ ^3/2 



1 ± 



1 x^/^ 



(8) 



(9) 



are the radii of the outer and inner Lindblad resonances for 
mode m. Waves are launched by the planet at the Lindblad 
resonances and propagate into r > r+ (with fc > 0) and into 
r < r_ (also with fc > 0). The sign of k is chosen such that 
the radial group velocity is directed away from the planet. 

^ Although torque calculations are sensitive to slight deviations 
from Keplerian rotation, e.g. resulting from radial pressure gra- 
dients, they are not important in the present analysis. 



2.3 Waves in the outer disc 

Consider the waves launched in the outer disc, i.e. the disc 
exterior to the planet's orbit. The precise phase of mode m 
for r > r+ is given by 



^m^^+Z A:(r')dr'-|-m((/>-t), 



(10) 



where the term 7r/4 is the phase shift associated with the 
resonance. This can be deduced from the behaviour of the 
Airy function, which describes the launched wave near the 
Lindblad resonance. 

We first obtain an estimate of the phase by a method 
equivalent to neglecting the term in the dispersion rela- 
tion. As noted above, this is appropriate when m is large. 
Approximating r± ~ 1, we have 



fc(r) dr ; 



-1 / 3/2 



1 dr. 



and therefore 

n 2m ( 



.3/2 



- In r — 1 
2 



1(0 -t). 



(11) 



(12) 



Constructive interference occurs near the curve <^ = V5(r, i) 
defined by 



t 



3/2 



3e V 
since we have 



■In r — 1 



$m ~ - + ra{4> - i^). 



(13) 



(14) 



We can, however, calculate the phase exactly by chang- 
ing to the variable x — r^^"^ and making use of the indefinite 
integral 



-(a + b) In [{x - a)'/' + {x - hf'^] - {abf^^ 

{a + h)x - 2ab - 2{aby/^{x - a^/^ix - bf^^ 



X In 



We write the exact solution as 



= - + m( 
4 



3e 



(15) 



(16) 



The last term measures the error in the approximate 
method, and will determine whether constructive interfer- 
ence really does occur on the curve 4> = V- The approximate 
method always slightly overestimates both the wavenumber 
k and the range of integration in the phase integral. It fol- 
lows that the residue Am has the properties 

dA„ 



A™ > 0, 



dr 



> 0. 



The expansion for large r is 
Am(r) = A„i(oo) - 



1 



2mr3/2 4mr3 



(17) 



(18) 



where (see Fig. 1) 

A„(oo) = m\n{2m) + {m'^ \n [m - (m^ - 1)^/^] .(19) 

A further property is that, for any fixed r. 
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Figure 1. The residue Am(oo). 



Figure 2. The residue Am in the outer disc. Modes m = 1, . . . , 10 
are shown from top to bottom. 



(1 



'3/2\ 



(20) 



(21) 



lim A„i = 0. 

m — 'oo 

In fact, for large m, 

Am ~ -^In \2e^^\ 
2m 

Waves of different m add constructively on the curve 
(j) — if provided that their phases lie within a range less 
than approximately vr. This certainly occurs for waves of suf- 
ficiently large m, because of property (po|). Constructive in- 
terference may fail for low values of m for which Am ^ 3ne/2, 
since such modes are out of phase with waves of high m. The 
properties ( |r^ ensure that the limit r — » cx) is the worst 
case. The residue is plotted in Fig. 2 for m = 1, . . . , 10. For 
e = 0.1, typical of protoplanetary discs, constructive inter- 
ference fails only for 1 < m < 2 in the limit of large r. 

Provided that a range of azimuthal wavenumbers is 
present, and there is no special selection rule (e.g. one that 
selects azimuthal wavenumbers that are all multiples of 2), 
there is no other curve on which constructive interference 
occurs consistently. The result is a one-armed spiral wake 
following the curve (f) = ip. 

The disturbance generated by a planet orbiting in the 
disc is dominated by azimuthal wavenumbers m ~ m* ~ 
l/(2e) (Goldreich & Tremaine 1980). Accordingly, the az- 
imuthal thickness of the wake is Sep « 2ii/m» « 47re. The 
radial thickness is smaller, Sr /r x e 6(f) x 47re^. 

2.4 Waves in the inner disc 

Waves launched in the inner disc behave very similarly. The 
precise phase of mode m for r < is given by 



$m = -+/ k{r')dr' + m{(t>-t), 



(22) 



where, again, > 0. We again write the exact solution in 
the form 



TT 2 

'^m = j+ m{(l) - f) + —Am{r) 



(23) 




Figure 3. The residue Am in the inner disc. Modes m = 1, . . . , 10 
are shown from top to bottom. 



where now 
2 



^-t+-[r^^-^-lnr-l 



In the inner disc the residue has the properties 
dA„ 



Am > 0, 



dr 



< 0. 



The limiting form for small r is 

Am(r) = -| [m - (m^ - 1)^^^] Inr + 0(1) 



(24) 



(25) 



(26) 



Constructive interference does eventually fail for all m in 
the limit r ^ 0. 
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The residue is plotted in Fig. 3 for m = 1, . . . , 10. The 
limiting form for large m is 

~ In [26^/^771 - l)] . (27) 



3 NUMERICAL CALCULATION 

To verify the hypothesis that the wake is formed through 
the constructive interference of Fourier modes, we have 
calculated the linear response of the disc numerically us- 
ing methods similar to Korycansky & Pollack (1993).q| We 
adopted a two-dimensional barotropic disc model equivalent 
to that described in Section and having a surface den- 
sity S oc r~^^^. Outgoing- wave boundary conditions were 
applied at r = 0.3 and r = 3. The potential of the planet 
was smoothed to simulate the non-zero vertical extent of 
the disc, with a smoothing length erp comparable to the 
semi-thickness H. The linear solutions for modes m = 1 to 
m — 100 were calculated and synthetized in real space. 

The results are shown in Figs 4 and 5, where we compare 
the predicted shape of the wake given by equations ( p^ 
and ( p^ ) with the numerical calculation. The agreement is 
excellent. The numerical calculation shows that the wake 
has some non-trivial internal structure, often consisting of 
both a trough and a larger peak. These details are likely to 
depend to some extent on the disc model. 

Although the waves are generated at Lindblad reso- 
nances, individual resonances cannot be observed in the 
complete solution because they overlap. The wake is gen- 
erated in a perfectly smooth manner. It is clearly estab- 
lished after the synthesis of the first 10 modes in the case 
of e = 0.05, or the first 5 modes in the case of e = 0.1. The 
addition of higher-order modes merely increases the finesse 
of the interference pattern. 



4 DISCUSSION 

We have considered the dynamical interaction between a 
planet of low mass and a two-dimensional gaseous disc in 
which it orbits. Although the planet generates disturbances 
of all azimuthal wavenumbers, the wave modes interfere con- 
structively on a unique curve and form a one-armed spi- 
ral wake. The waveform based on linear theory (see Figs 4 
and 5) is in good agreement with that found in non-linear 
planet-disc simulations (e.g. Artymowicz 2001). 

The formation of a coherent structure by constructive 
interference is reminiscent of the Kelvin wedge produced in 
the wake of a ship. In the present case the effect is more 
complete, as almost all modes interfere constructively, rather 
than a band of wavenumbers. 

Recently, Goodman & Rafikov (2001) calculated the lin- 
ear wake produced by a terrestrial planet in a local, shearing- 
sheet model of a two-dimensional disc. As they noted, the 
formation of the wake enhances the amplitude of the distur- 
bance in a thin structure and therefore increases the likeli- 
hood of non-linear effects such as shock formation. 

^ We caution the reader of that paper of a number of errors, 
especially on p. 164. 



The precise relation between a linear wake and the spi- 
ral shocks observed in numerical simulations involving plan- 
ets of larger mass (e.g. Lubow et al. 1999) is not entirely 
clear, however. In such simulations a dominant one-armed 
shock is formed, which follows a similar curve to the linear 
prediction, but other features are also present. The flow near 
the planet is also different in the non-linear regime. 

Although most treatments of the planet-disc interac- 
tion have used a two-dimensional description of the disc, 
this is in fact difficult to justify. First, the wave modes in a 
three-dimensional disc are, in general, different from those 
in a two-dimensional disc. If the disc is vertically isothermal, 
there does exist a two-dimensional mode with the same dis- 
persion relation as equation ^ (Lubow & Pringle 1993) . In 
that wake would be formed in much the same way. 

In a vertically thermally stratified disc, however, the waves 
generated by the planet have a quite different dispersion re- 
lation. For example, in a polytropic model, which represents 
a highly optically thick disc with vertically distributed en- 
ergy dissipation, the equivalent mode behaves like a surface 
gravity wave having an approximate dispersion relation 

[m(n - Qp)f « gk (28) 

sufficiently far from the planet, where g = Q^H is the ver- 
tical gravitational acceleration at the surface of the disc 
(Ogilvie 1998; Lubow & Ogilvie 1998). Spiral waves with 
different values of m then have different pitch angles and 
there is no possibility of consistent constructive interference 
on a curve.n We therefore anticipate that the disturbance 
generated by a planet in a thermally stratified disc will have 
a quite different structure. 

In addition, the flow near the planet is likely to be three- 
dimensional if the radius of the planet's Roche lobe is com- 
parable to or less than the semi-thickness of the disc. There- 
fore, although this analysis may explain some features ob- 
served in numerical simulations, the realities of the planet- 
disc interaction are likely to be more complicated. 

We remark that an accurate linear calculation of the 
planetary wake, and of the associated migration rate, in a 
thermally stratifled, three-dimensional disc would be a de- 
manding numerical problem, and has not been carried out. 

In summary, the one-armed wake is a consequence of 
having a spectrum of two-dimensional waves with differ- 
ent azimuthal wavenumbers m and approximately matching 
phases. Any mechanism capable of producing such waves 
will result in a similar one-armed wake. We have shown that 
resonantly launched waves in a two-dimensional gaseous disc 
naturally meet this criterion. 
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in a ship wake, the analogy with the Kelvin wedge breaks down 
when the shear in the disc is taken into account. 
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Figure 4. Left: Predicted shape of tlie spiral wake for e = 0.05, based on equations ( |l3| ) and (p4|). Tlie dotted line represents the 
corotation circle, r = 1. The planet is located at (1,0), and the outer radius plotted is r = 3. Right: Numerically calculated spiral wake 
for € = 0.05. The enthalpy perturbation is plotted using a linear grey-scale from negative (black) to positive (white). The maximum 
intensity corresponds to a fractional surface density perturbation, at r = 1, of 10*(Mp/M). 
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